Superfiuid equation of state of cold fermionic gases in the Bose-Einstein regime 
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We present an exact many-body theory of ultracold fermionic gases for the Bose-Einstein conden- 
sation (BEC) regime of the BEC-BCS crossover. This is a purely fermionic approach which treats 
explicitely and systematically the dimers formed in the BEC regime as made of two fermions. We 
consider specifically the zero temperature case and calculate the first terms of the expansion of the 
, chemical potential in powers of the density n. We derive first the mean-field contribution, which has 

the expected standard expression when it is written in terms of the dimer-dimer scattering length 
a_M- We go next in the expansion to the Lee- Huang- Yang order, proportional to n'^^'^ . We find the 
far less obvious result that it retains also the same expression in terms of am as for elementary 
bosons. The composite nature of the dimers appears only in the next term proportional to n^. 

I. INTRODUCTION 

^ ' The recent experimental developments in ultracold Fermi gases and in particular the experimental realization of 
the BEC-BCS crossover through Feshbach resonance [l[ have arised in a field where theoretical investigations started 
a long time ago. Indeed fairly soon after the publication of the Bardeen-Cooper-Shrieffer (BCS) theory 0, the 
extension to the neutral Fermi liquid ■^He was considered 3. Although the case of liquid ''He is more complex than 
the one of ultracold gases, it is nevertheless analogous to the situation found on the BCS side of the crossover, in 
^ ■ the strong interaction regime. Similarly the possibility of Bose-Einstein condensation (BEC) of composite bosons was 
considered quite early for excitons in semiconductors [1] , where the composite nature of these bosons is expected to 
be an essential physical feature. This is the situation met in the BEC range of the BEC-BCS crossover. 

Deep theoretical investigation of this last situation was made not long afterwards by Keldysh and Kozlov 0], 
'"^ ■ following a seminal work by Popov Q which considered a physical model with short range potentials, much closer to 
^ I the ultracold gas case. In both cases, due to the formal similarity between the Bogoliubov and the BCS frameworks, 
O . the theoretical treatment relied on the fact that in the dilute limit the formalism leads directly to the Schrodinger 
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equation for the molecular state corresponding to the coinposite boson. This provided an anticipation of the study 
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of the whole BEC-BCS crossover. Independently Eagles Q, in the course of a study of superconductivity in doped 
semiconductors (where the electron gas has a low concentration compared to standard metals), was led to investigate 
^ ' the extrapolation of the simple BCS formalism toward the dilute regime where pair formation would take place, due 
to the attractive interaction, and Bose condensation of these pairs should occur at lower temperatures. 

In the context of cold gases, the physics of the crossover was considered by Leggett Q, who stressed the physical 
interest of considering Cooper pairs as giant molecules, emphasized the change of physical regime when the chemical 
potential goes through zero, and introduced the scattering length in the BCS formulation. The matter was taken 
up by Nozieres and Schmitt-Rink Q who studied the smooth evolution of the critical temperature in the crossover 
within the simplest T— matrix approximation, showing in particular how the critical temperature of the ideal Bose 
[ gas is recovered in the strong coupling limit. Sa de Melo, Randeria and Engelbrecht [10] addressed the crossover with 
the functional integral formalism. 
. ^ ^ As it has often been stressed, the physical situation found in this BEC-BCS crossover is extremely interesting 
. since (at least for wide Feshbach resonances (T], [ll|) the scattering length a of the different fermions with mass m is 
^ ■ enough to fully describe their interaction, while one has still to deal with the full complexity of a strongly interacting 
system. As a result the theoretical problem has not been solved exactly and one has to rely on various approximate 
schemes the simplest one being the straight BCS formalism (often called BCS-MF) which is known to give a 
correct description at zero temperature in both the extreme BEC and BCS limits a — > 0±. The imperfections of the 
approximate treatments are clearer in the BEC regime. Indeed if one wants to describe properly the departure from 
the ideal Bose gas of dimers, one has to provide a correct description of the dimer-dimer interaction, in order to find 
the expected physics of a weakly interacting Bose gas. As stressed by Keldysh and Kozlov [s^ the composite nature 
of the dimers has already to be properly taken into account at this stage in order to obtain correct results. 

A first ste p in this direction for the specific case of short range interactions, relevant for cold gases, was made by 
Haussmann who obtained, both in the normal and in the superfiuid state (where this can be found naturally 
out of the BCS ansatz), that the dimer-dimer scattering length om is related to the fermion scattering length a by 
um = 2a. This result was also obtained in Ref.jT^ in the superfiuid state by other methods. The level of approximation 
corresponding to this result turns out to be the Born approximation for the dimer-dimer scattering. At the same level 
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of approximation Fieri and Strinati [13| derived quite recently the Gross-Pitaevskii equation from the Bogohubov- 
de Gennes equations. However this level of approximation was markedly improved by Fieri and Strinati [14j who 
performed the calculation at the level of a T— matrix approximation and found om — 0.75 a. Nevertheless their 
calculation does not include all the processes resulting from the existence of the fermions making up these composite 
bosons. The exact numerical result was obtained recently by Petrov, Salomon and Shlyapnikov [151] who basically 
solved the relevant Schrodinger equation for the four-body problem and found for this dimer-dimer scattering length 
qm — 0.60 a. Thereafter it was shown fl6] how this same problem could be exactly formulated and solved within 
the methods of quantum field theory, with naturally exactly the same result. This last work is quite relevant to the 
present paper, and represents the basis out of which we will work. This is naturally linked to the fact that field 
theoretic methods are a convenient, and probably necessary, framework for an exact formulation and solution of this 
superfluid strong coupling problem raised by the BEC-BCS crossover. 

In this paper, we provide an exact approach to the BEC-BCS crossover problem from the BEC side by presenting 
an exact fermionic theory of a BEC superfluid of composite bosons in the low density range. The two interacting 
fermion species we deal with are actually the two lowest energy hyperfine states used experimentally [l| in ultracold 
gases of ^Li or ^'^K, and as it is usually done we will for convenience refer to them as "spin up" and "spin down" 
states. Specifically working in the low density range means that we proceed basically to an expansion in powers of the 
gas density. Otherwise our framework is completely general. In the present paper we deal with the lower orders in 
this expansion, but nothing in principle prevents our approach to be extended to higher orders, although admittedly 
this will require much more work. Similarly we will restrict ourselves to the T — situation, which brings strong 
simplifications as we will see because in this case there are no fermions at all in the normal state. However it is 
quite possible to extend our approach to non zero temperature. Finally we deal here with thermodynamics, but our 
framework allows us naturally to calculate dynamical quantities. 

As we mentionned our approach is in principle a low density expansion. However the density n is not a convenient 
basic parameter since it rather appears in the formalism as the result of a calculation. On the other hand it is known 
(and we will see it explicitely below) that in the simple BCS approximation for the BEC limit, which as mentionned 
above corresponds to a Born approximation at the level of the mean field term, the density n is proportional to A^, 
where A is the gap parameter of the BCS theory. Hence a convenient way to proceed effectively to a low density 
expansion is rather to expand in powers of A. However this quantity is specific of the BCS approximation, which has 
no general validity. But it is directly related to the anomalous self-energy A(p) (where p is an energy- momentum four- 
vector) which is a completely general microscopic quantity characteristic of a superfluid system. Hence our method 
will specifically be an expansion of the general equations in powers of the anomalous self-energy A{p). 

Here we will calculate the equation of state, or more precisely the dependence of the chemical potential /i on the 
density n (a short account of our calculation has already been published [l3|)- Actually this is rather the reciprocal of 
this function which comes naturally out of the calculation, since /i enters the Green's functions while n is obtained at 
the end of the calculation, after elimination of A(p) which comes in the intermediate steps. The lowest order term in 
the expression of /i, beyond the trivial term of half the molecular binding energy — l/(2ma^), is the mean-field term. 
One has to note that, even at this stage, the composite nature of the bosons we deal with is entering as noted by 
Keldysh and Kozlov [5] . However in our case it is lumped into the full dimer-dimer scattering length om ■ This result 
is physically natural since, beyond the zeroth order term — l/(2TOa'^) in the chemical potential [1^ corresponding to 
an isolated molecule, we expect in a density expansion to find the physics of two coexisting molecules, which should be 
fully describable by gm for its static properties, corresponding to the zero energy scattering properties. Nevertheless, 
as far as we know, this quite reasonable result has always been in the literature more or less taken for granted or 
cursorily obtained, but not fully explicitely derived from a microscopic theory for composite bosons. Hence our first 
step will be to obtain this fully explicit derivation. 

However we are naturally essentially interested by the next order term, at which we will also limit our present 
calculation. From the theory of weakly interacting elementary bosons this term is proportional to n'^/^, as first 
obtained by Lee and Yang [19|. We refer to this term as the LHY term. This term has then been rederived by 
Beliaev [20] making use of field theoretic methods. Naturally it is reasonable to expect that a similar term arises for 
composite bosons. What is however less clear is that the coefficient is just the same as for elementary bosons, i.e. that 
it can be simply expressed in the same way in terms of gm- Indeed it is clear that, at some stage in the expansion, 
the composite nature of the bosons will enter by other quantities than those describing the physics of elementary 
bosons. This is quite obvious physically since, by going to higher densities (which corresponds to take into account 
higher order terms in the expansion), we will go toward unitarity and find physical properties linked to the existence 
of Fermi seas. Obviously this can not be obtained from the physics of elementary bosons. Hence it is conceivable 
that additional processes contribute to the n^^^ term for composite bosons. Nevertheless we will find at the end of 
our paper that the nP^^ term is indeed identical to the LHY term provided the dimer-dimer scattering length is used. 
The fact that this is not obvious is immediately realized from the number and the complexity of our steps, which do 
not map systematically on the elementary boson derivation. 
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On the other hand this result makes sense physically when it is realized that the LHY term is directly linked to the 
existence of the gapless collective mode. A standard perturbation expansion should produce a rt^ term as the next 
term after the mean-field one. The existence of the LHY term is linked to the presence of these gapless elementary 
excitations which produce a singularity in the expansion, leading formally to a divergence in the coefficient of the 
term. Obtaining the proper result is equivalent to the resummation of a partial series to all order in perturbations in 
order to get rid of the singularity. Since these steps involve only low energies, they do not test the internal structure of 
the bosons which comes into play when energies of order of the binding energy 1/ma? are involved. In such a case it is 
reasonable that no difference appears between elementary and composite bosons. On the other hand the calculation 
of the coefficient of the regular n? term in the expansion involves the consideration of energies of order of the binding 
energy, as we will see explicitely, making explicit the composite nature of the bosons. This implies that at the level 
of this term one finds a result different from the case of elementary bosons. Note that this result has been previously 
assumed to be correct in calculations of collective mode frequencies [2l|, . On one hand this was in agreement with 
Monte-Carlo calculations [l^, and on the other hand this has been supported by very recent experiments [23|. 

At the level of handling the collective mode, very strong analogies with a purely bosonic approach, in particular 
with Bcliaev's work [205 ^'"^ appear for obvious physical reasons, although it can be seen that our treatment bypasses 
a good deal of Beliaev's detailed diagrammatic analysis. We will not try to put these analogies away, but rather 
make use of these similarities and correspondances to make easier the physical and technical understanding of our 
calculation. However we will never make an actual use of a bosonic approach and all our theory is a purely fermionic 
theory, the essential physical feature being that we deal with fermion physics but with negative chemical potential 
/i < 0. In particular we often use the word dimer to describe two fermions with spin up and spin down. But since the 
theory contains, in the corresponding propagator, not only the bound state of these fermions but also the scattering 
states, it should be taken as a simple convenient wording, although it is naturally physically suggestive. More generally 
and formally the anomalous self-energies A(p) and A*(p) play clearly a role analogous to (6o) and {b\) for bosonic 
superfiuidity. 

Our paper is organized according to the above introduction. After introducing briefiy the general formalism in 
the next section [Hi we proceed to the calculation of the normal and anomalous self-energies in section IIIII to lowest 
order in A(p), which leads us to the mean field term in the expression of the chemical potential. Then in section HVl 
we proceed to discuss the collective mode which plays an essential role in the calculation of the LHY term. Then 
the calculation of the contribution of this collective mode to the normal and anomalous self-energies is performed in 
section |Vl at the end of which the LHY contribution to the chemical potential is obtained. A number of technical 
details are treated in appendices, to try to make more palatable the main part of this paper which is already fairly 
heavy. 

II. GENERAL 

Our work deals only with the T = situation. Hence it would be most natural to use the standard T = Green's 
function formalism and naturally all the content of our paper can be reproduced with this formalism. However 
one has in this case to be careful in dealing with real frequencies, and specify precisely the location, with respect 
to the real axis, of singularities in the corresponding complex plane. This is often left implicit. A simple way to 
clarify this question is to deal with complex frequencies. But this is actually what one has to do naturally when 
one works at T 7^ 0, where the frequencies run on the imaginary frequency axis. In the T — s- limit the discrete 
sums over Matsubara frequencies turn into integrals over frequencies running on the imaginary axis (the link with the 
T = formalism is merely obtained by continuing analytically the frequencies toward the real frequency axis, that is 
performing a so-called Wick rotation on the frequency). This leads to somewhat clearer results. Hence, despite its 
somewhat artificial character, we will work quite specifically with this T — Q limit of the T formalism. Actually 
this will not appear in many our expressions because of their formal nature and for this reason it does not make 
any real problem. The essential point is that, in this way, we will always have perfectly regular expressions, with an 
effective handling which is quite clear and easy. 

We deal only with the formalism of fermionic superfluidity, as it is known in superconductivity 25] and superfiuid 
■^Hc 26]. This entails the use of both the 'normal' fermionic Green's functions G'||(p) — Gii{p) = G{p), where we 
make use of the four- vector notation p = {p, v}, and their 'anomalous' counterpart F{p) and F~^{p), where we make 
use of standard notations [1^. This could he written in a compact Nambu formalism, but for clarity we prefer to write 
the normal and anomalous parts explicitely, with slightly different sign conventions. In practice we will always consider 
in the following, when we calculate G(p), that we deal with G||(p), without writing it explicitely for simplicity. The 
existence (i.e. the fact that they are non zero) of these anomalous Green's function is the basic ingredient of fermionic 
superfluidity in our formalism. This is quite clear physically since they are related to the condensate. 



FIG. 1: Diagrammatic representation of Dyson's equations Ea. (lll2|l 



We start by writing Dyson's 




(1) 
(2) 



their corresponding diagrammatic representation being given in Fig. [TJ Here Gq(p) = [iv — ■^p]^^ is the free atom 
Green's function, with — — fJ- the atom kinetic energy ep — p^/2m measured from the chemical potential /i. 
The quantities and A(p) are respectively the normal and the anomalous self-energies. By definition is the 
sum of all diagrams, with one entering p line and one outgoing p line, which can not separated into two disconnected 
parts by cutting a single Go{p) line (or Go{—p) line). A{p) has the same definition except that it has two entering 
lines, one (p, t) line and one {—p, I) line. We note that A(p) depends naturally in the general case on momentum and 
frequency. In the second equation, for which begins with an outgoing (— p, J,) line and finishes with an outgoing 

{p, t) line, it would look more symmetrical to introduce a corresponding anomalous self-energy A'^{p), analogous to 
A(p), but with two outgoing lines, one {—p, I) line and one (j>, t) line. However it can be related directly to A*(p) by 
time reversal, as we have written explicitely. All formulae in the following turn out to be simpler to interpret when 
this substitution is done, and we will follow this standard habit. We note also that, since we deal with an equilibrium 
situation, where the gas is homogeneous, we do not need to introduce Green's functions with different entering and 
outgoing momenta, as it would be required in the case of an inhomogeneous system. 

Dyson's equations Eq. ([T]) can be understood as a way of ordering the full perturbation expansion, by gathering 
terms into and A(p). However it is interesting to note that this can be done in slightly different ways. For 

example we can gather in Eq. (jT]) all the terms where only the free Green's function Go{p) and the normal self- 
energy appear, that is the terms where no anomalous terms are present. This leads to introduce Go{p) — 



Go{p) + Go{p)nP)Go(.p) + Go(p)S(p)Go(p)S(p)Go(p) + that is g^\p) = G^\p) ~ In this way, by making 



use of Goip) instead of Go{p), we are only left to write in Dyson's equations the anomalous terms, that is to rewrite 
these equations as: 



It is indeed checked easily that the algebraic solution of Eq. (Ill2p . as well as the one of Eq. ([Sj H]), is indeed: 



This result could have been obtained directly in the following way. In the definition we have taken above for the 
self-energy, we had excluded diagrams made of two blocks linked by a single propagator Gq{p) (corresponding to an 
atom propagating to the right, as for the definition of the normal self-energy) or Gq{—p) (corresponding to an atom 
propagating to the left, as it is produced by the anomalous self-energies). However we could take another definition 
where we would exclude only blocks linked by a rightward propagator, just as in the normal state. In this case the 
formal expression of G{p) in terms of the self-energy, i.e. Gq^{p) = Gq^{p) — 5(p), is just the same as in the normal 
state, as it is indeed the case in Eq. But the expression of the self-energy has been naturally modified from S(p) 
to S{p). Indeed we have to include now terms containing leftward propagators Go{—p). But it is easily realized that 
these new terms contain at most a product A(p)A*(p), corresponding to the presence of anomalous self-energies of 
opposite types at the extremities of the block. These two anomalous self-energies can be separated by any number of 
leftward propagators Go{—p), with normal self-energies E(— p). The sum of all these produces a factor Go{—p), which 
leads to the relation: 



G{p)^go{p) + go{p)A{p)F+{p) 

F+{p) = -go{-p)A*{p)G{p) 



(3) 
(4) 




(5) 



(6) 



s{p)^j:{p)^A{p)g,{~p)A*{p) 



(7) 



in agreeement with Eq. 
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FIG. 2: Diagrammatic representation of 5'-^-'(p) 

III. A EXPANSION: LOWEST ORDER 

While they are completely general, the considerations presented in the preceding section are not of much practical 
use since we have no expressions for the normal and anomalous self-energies. In order to obtain explicit expressions, 
we will now begin to proceed along the lines outlined in the introduction, namely to perform a systematic expansion 
in powers of A. In the present section, as an introduction and as a first step, we will only consider the lowest order 
in this expansion. The next order will then be addressed in the next section. In the first subsection we will consider 
the expansion of the normal self-energy, which turns out to be the easier one. Then the anomalous self-energy will be 
handled along the same lines. 

A. Normal self-energy 

As we have indicated in the introduction, and as wc will check explicitely below, the expansion in powers of A has 
the character of an expansion in powers of the density n since we will find n ~ A^ . Hence the zeroth order in powers 
of A corresponds to an infinitely dilute gas for which the normal self-energy is naturally zero. Since by its definition 
this normal self-energy contains only processes which conserve globally the number of particles, the lowest order term 
Yi^'^\p) in the normal self-energy must contain both A and A*, because A implies from its definition the annihilation 
of two atoms (going physically in the condensate), while A* corresponds to the creation of two atoms (coming out 
of the condensate) . Hence this lowest order term in the normal self-energy is the sum of all diagrams which contain 
A(p') and A*(p"), and which have one entering line p and one outgoing line p, the rest of the diagrams containing 
only (since we proceed to an expansion) A = quantities, that is normal state quantities (see Fig. [5]). We call the 
sum of all these normal state diagrams the 'normal part' of S(^)(p). This is just I](^^(p), except for the A{p') and 
A*{p") terms. 

Explicitely this means that this 'normal part' of the diagrams contains only any number of bare propagators Go(p) 
and of interactions lines. Moreover the two atom lines coming out of A*(p")(with opposite spins) are not allowed to 
interact immediately after entering the 'normal part' of the diagram, since such an interaction is already included in 
A*(j>") (as we will see explicitely below) and allowing it would amount to double counting of diagrams. Similarly 
the two atom lines entering A(p') are not allowed to interact immediately before. As it happens, such a normal state 
set of diagrams has already been essentially considered in Ref. JjjJ. Indeed in this paper the scattering of a single 
atom by a dimer (made of atoms with opposite spins) has been calculated. The fact that a dimer is entering the 
diagram implied that the two corresponding atom lines can not interact immediately in the diagram, since such an 
interaction is already taken into account in the dimer propagator and this would be again double counting. There 
is a similar requirement for the two lines corresponding to the outgoing dimer. We see that the requirements on the 
normal part of E'-^-'(p) are just the same as in the above paper. In this Ref.[l6j the sum of all these diagrams was 
denoted T3{pi,p2; P) where pi (respectively P2) is the four- vector of the entering (respectively outgoing) atom, and 
P is the total four- vector of the single atom and the dimer. In our statement we apparently overlook the fact that 
T3{pi,p2] P) contains only three atom lines, which go from the beginning to the end of the diagram, with various 
interactions between them in the diagram, since we deal with the vacuum scattering of an atom and a dimer. In 
contrast we could think that, in addition to these lines, the 'normal part' of I](^^(p) may contain any number of loops 
of atom lines. However it is clear physically that these loops are not allowed. Indeed if we think of this 'normal 
part' as representing, in a time representation, a succession of processes due to interaction, the appearance of a loop 
corresponds to the creation of a particle-hole pair. While such processes can always occur at non zero temperature, 
in our T = situation they are possible only if the chemical potential /i is positive, thereby allowing the existence 
of a Fermi sea in which the hole can be created. But in our dilute regime we have fj, < and these processes are 
forbidden. This argument is developped more technically in Appendix [BJ Hence it is correct to identify the 'normal 
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part' of E(2)(p) with since both contain only three atom hnes, subject to the same restrictions indicated above. 

Nevertheless there is still a small difference in the general case between the 'normal part' we want and T^. In the 
definition of T3 one has to sum over all possible four-vectors for the two lines making up the entering dimer, with 
the only restriction that their sum is fixed at P — pi . A similar summation applies for the two lines making up the 
outgoing dimcr. On the other hand in our 'normal part' the two lines (p',T) a-nd (-p',i) entering in A{p') have a total 
zero four-vector, and because A(p') depends on p', when we sum over p' , there is a weight A(p') which depends on 
p', in contrast with the case of T3. Nevertheless, as shown in Appendix [XI it is possible to express in the general case 
our self-energy in terms of T3. But in our specific problem, this is an unnecessary complication since we will see that, 
at the order of the expansion we are considering, the dependence of A(p') on p' can be neglected, and A(p') can be 
considered as a constant which we denote merely as A. Hence we can factor A out. Similar considerations apply to 
A*(p") and the entering dimer lines. 

There is still a very important point to consider before writing our expression for S(^)(p). Among all the diagrams 
we have considered contributing to S(^)(p), there is necessarily a diagram which gives a contribution — AGo(— p)A* 
because it satisfies all our requirements. Indeed this is just what is coming from the lowest order diagram in T3, which 
gives a contribution — G'o(— p) to T3 and corresponds merely to the Born approximation. This diagram is necessarily 
present because we considered it already in the preceding section [III when we were dealing with the general expansion 
of G{p) in powers of S(p) and A(p). This is just the last term in Eq.([6]), when we replace Goi—p) by Go(— p) (since we 
want only 'normal state' diagrams) and take A(p) independent of p. However, since it is already explicitely present in 
Eq.®, this implies automatically that it is reducible (and indeed it is obviously reducible) and for this reason should 
not be included in S(^)(p) (otherwise we would double count this diagram). On the other hand this is clearly the 
only reducible contribution coming from T3, as seen directly from the diagrammatic expansion [l6j of T3, or because 
any such reducible contribution has to appear explicitely from Eq.([S]) as explained above. Accordingly we have to 
subtract this contribution out and we obtain finally: 

S(2)(p) = |A|2T3(p,p;p) + \AfGo{-p) (8) 

We can also rewrite this result in terms of our self-energy S as: 

5(^)(p) = |A|2T3(p,p;p) (9) 

This result is also directly quite natural since, when we evaluate S, the diagrams with Go{—p) are not considered as 
reducible, as we have seen, and accordingly — AGo(— p)A* should not be subtracted. 

We are now in position to calculate the single spin atomic density n, which is given quite generally by: 

7^Pf^G(p)e-^ (10) 

If G is replaced by Gq, one checks easily that n — 0, which is physically obvious since in this case we deal with a non 
interacting fermion gas, at T = 0, and chemical potential /i < 0. There is no Fermi sea (which arises only for fj, > 0) 
and there are no fermions. Hence we can replace completely generally G(p) by G(p) — Go(p) = Go(p)iS(p)G(p) in 
Eq. pUj) . Since we work in this section at the lowest order in A, which means at order A^ in the present case, we can 
replace S{p) by iS^^^(p), and G(p) by Go(p). This leads to: 

"HA|»/|Er£°°|:;^3(P.P.rt|Grf (11) 

where, as in Eq. pI7|) . p is for {p, i^}. 

At this stage it seems that, in order to calculate the density, we are required to have a deep knowledge of T3. This 
seems physically reasonable since, in Eq. (|lip . T3 describes the scattering of an atom with the condensate which has 
clearly to be taken into account. However the actual situation turns out to be much simpler. In order to calculate 
the integral over v in Eq. (jlip . it is convenient to close the contour in the upper complex half- plane Imz/ > (which 
corresponds to Re (w) < 0), since the pole of Go(p) = [iv — £,p]~^ is found for Imi' < (since = ep -I- |^| > 0). 
However it is shown in Appendix [Cl that, except for the Born contribution, T3(p,p;p) is analytical in the upper 
complex half-plane and does not contribute to the integral. The result comes entirely from the Born contribution. 
This means that surprisingly the contribution of I]''^-'(p) to n is zero, and that the only non zero contribution comes 
from the explicit BCS term — |ApGo(— p) in Eq. ([7]). In other words the approximate BCS formalism gives the exact 
result at the lowest order in our A expansion. It reads: 

-1^1^ /i?^/ri^ 
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FIG. 3: The two sets of diagrams making up the anomalous self-energy. By definition the second set excludes the diagrams which 
are already included in the first set. In this second set the full propagators drawn may be normal or anomalous (accordingly we 
have not drawn the second arrow on these propagators). Similarly the vertex corresponding to the rectangular box has normal 
and anomalous parts. 



If we make use m this result of the zeroth order value of the fcrmion chemical potential /i — —Ei,/2 = —\/2ma?', we 
end up with: 

™2 „ I A 12 

(13) 



Stt 

Actually this relation |Ap = (87r/m^a) n is precisely obtained from the dimer propagator T2(P) (see section ITV]) when 
it is evaluated near its pole and with the purely bosonic propagator factor merely replaced by n. 



B. Anomalous self-energy 



We proceed now in a somewhat similar way for the anomalous self-energy. The resulting equation is the general- 
ization of the standard 'gap equation' of the conventional BCS theory and it has the same character of being finally a 
self-consistent equation for A(p). In writing an expression for A(p), we can in full generality separate all contributing 
diagrams in two sets. In the first set, for both of the two entering fermionic lines, the first interaction is between these 
lines themselves. The second set corresponds to the opposite case, where for at least one of these two entering lines 
the first interaction is with a line which is not the other entering fermionic line (see Figl3]). This separation of the 
diagrams in two classes is completely analogous to the separation of the diagrams, coming in the standard self-energy 
(25} . into an Hartree-Fock like set and another set involving the full vertex. 

The set of diagrams where the two fermionic lines interact can be immediately written in terms of F{p), which is 
analogous (and actually directly related) to F^{p), except that it has, instead of outgoing lines, an ingoing {p, |) line 
and an ingoing {—p, i) line. It satisfies a Dyson's equation analogous to Eq. pi2p . namely: 

G{p) = Goip) + G{pMp)Go{p) - F{p)A*{p)Go{p) (14) 
F{p) = F{p)i:{-p)Go{-p) + G{p)A{p)Go{-p) (15) 

Introducing the Fourier transform V{q) of the bare instantaneous interaction between |-atom and |-atom, this first 
contribution to the anomalous self-energy is: 

6iip) = - J j^Vip-k)Fik,u;) (16) 

which has naturally just the form of the standard BCS term. 

Let us now switch to the other set of diagrams, for which we proceed to an expansion similar to the one carried 
out in the preceding subsection. Because of particle conservation this expansion involves now odd powers in (A, A*). 
We consider first the term proportional to A(p'). The 'normal part' of this diagram contains only one t line and 
one I line (with any number of interactions), since additional lines are not allowed as discussed above. Moreover, 
because the interaction is instantaneous, there is no possibility of crossing of interactions when the diagram is drawn 
in time representation. Finally our lines are not allowed to go backward in time since this would correspond to the 
creation of a particle-hole pair, which is again not allowed in our case where we have a negative chemical potential 
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at T = (by contrast, for /i > 0, such processes could naturally occur). In summary the only possible diagrams 
are those of repeated interactions between the t and the | atom. This corresponds merely to the propagation of a 
normal dimer. However these diagrams are already generated by the iteration of the term (fT6|) when F{k) is written 
as Ga{p)A{p)Ga{—p) from Eq. ^TS]) at lowest order, and they should not be double-counted. Hence there is no term 
of order A in this set, and our expansion starts with terms containing A{p')A*{p")A{p"'). 

The analysis is now somewhat similar to the one made in the preceding subsection IHIAI To obtain the above third 
order term in the equation for A{p), we have to write, together with the factor A(p')A*{p")A{p"'), all the 'normal 
state' diagrams which have, in addition to the two entering lines (p, t) and (— p, i), another couple of entering lines 
coming from A*{p"). Moreover, we have two couples of outgoing lines, going into A(p') and A{p"'). Just as in the 
preceding subsection, we have a 'no first interaction' restriction for all these couples of lines, including the (p, "f ) and 
{—p, I) ones. As discussed above, in addition to these ingoing and outgoing ones with their continuations, no other 
lines can appear in these diagrams since they form loops, which are forbidden since we are at T = and /i < 0. We 
find that such normal state diagrams have already been essentially encountered in Ref. [l6j , where their sum has been 
called (72;P2, -P)- Here qi and q2 are for the entering fermionic lines (corresponding to p and — p in our case), 

the total four- vector for the ingoing and outgoing lines being 2P. Finally there is, in the definition of $, one entering 
2P — qi — q2 dimer, while there are two outgoing dimers P + P2 and P — p2- 

Just as in the preceding section, we should take into account that, in the general case, we have the weight 
A{p')A*{p")A{p"') which depends on the four- vectors of the atom propagators making up the dimers, while in 
the definition of <i> the summation over these variables is already performed, with only their sums 2P — qi — q2, P +P2 
and P — P2 being fixed. This problem could nevertheless be circumvented and an expression for the 'normal part' 
obtained in terms of $, in a way analogous to Appendix]^ However, just as in subsection IIII Al this is not necessary. 
Indeed from the lowest order term Si{p) given by Eg. p^ . we see that A{p) does not depend at lowest order on the 
frequency variable v. Moreover, since we deal with a bare interaction potential V^(r) which has a very short range, 
its Fourier transform l^(q) is a constant, except for very large wave vectors. But, as we will see, these wave vectors 
will not come into play because all the integrals we deal with converge at a much shorter wave vector. Hence we find 
that, at lowest order, A(p) turns out to be independent of p, just as in standard BCS theory. Now, when we want 
to calculate the third order term, we may insert in it the expression of A(p) to lowest order, which is just a constant 
which we denote again as A. Accordingly the above problem disappears and the 'normal part' of our diagrams can 
be expressed directly in terms of ^{qi,q2;p2, P)- 

Nevertheless we have to take into account, just as in subsection IIII Al that some of the diagrams appearing in 
A |Ap$ are actually not irreducible, and should not be included in our expression for A{p). Hence we have to 
subtract them out of A |Ap$. We denote $' what remains of $ after this subtraction. This will be taken care of just 
below. Gathering our results, we end up with the following (self-consistent) equation: 

A(p) = - 1 ^ V{p - k) F(k, tu;) + ^A I Ap $'(p, -p; 0, 0) (17) 

The factor 1/2 in front of (f>' is a general topological factor (it is present even if we do not take A as constant). Indeed 
exchanging A(p') and A(p"') corresponds to a mere change of variable and not to a different diagram. On the other 
hand, by making use of $'((71, q2',P2,P), any given such diagram would appear twice (even if $ is actually unchanged 
when we exchange the bosonic variables of the two outgoing dimers, which amounts to change P2 into — P2). Hence 
the factor 1/2 is required to avoid this double counting. We note that clearly A(p) depends indeed on p in general, 
since this is already explicitely the case in Eq.(fT7|) where the third order term depends on p through $'(p, — p; 0, 0). 

Let us now be specific about the difference between $ and due to the reducible diagrams contained in A |Ap$. 
One can obtain them by looking at the diagrammatic expansion [iBl of $. However it is easier to notice, as indicated 
above, that they appear automatically explicitely from Ea. (ll4ll5p . when their solution are expanded in order to 
obtain the complete perturbative series expansion of G'(p) and F{p) in powers of S and A, because this expansion 
provides precisely all the reducible contributions in terms of the irreducible self-energies S and A. When we expand 
the solution for F{p): 

F{P) = r^^^^^ (18) 

[Go-i(p) - E(p)][Go-i(-p) - E(-p)] + A(p)A*(p) 

to third order, we find: 

F(p) = Go(p) \a{p) + E(2)(p)Go(p)A(p) + Aip)Go{~p)^^^H-p) - A(p)Go(-p)A*(p)Go(p)A(p)l Go(-p) (19) 



where for consistency we have to take only the second order contribution S'^^^(p) (given by Eq.®) to E(p). The third 
order terms in the bracket are clearly the reducible contributions to A | Ap"I>(p, — p; 0, 0), as it can be seen directly. 




T 




+ 






FIG. 4: Reducible contributions to the anomalous self-energy from A |A|^ ^(p, — p; 0, 0). In the first two terms is obtained 
from T3 by excluding the Born term, since it is taken into account explicitely by the third term. 

They are displayed diagrammatically in Fig[4l In particular the last term of the bracket comes from the Born term 
of —p; 0, 0). Hence we have explicitely: 



-A|A|2cO'(p,_p;0,0) 



-A|A|2$(p, -p;0,0) 

I](2)(p)Go(p)A(p) + A{p)Go{-p)^^'H-p) - Aip)Go{~p)A*ip)Goip)A{p) 



(20) 



However it is now more convenient to work with <1> rather than Hence we add to both members of Eq.([T7]) the 
reducible contributions. After multiphcation by Go{p)Go{—p) and taking into account Eq. p9|) . we end up with: 



F{p) = Go(p)Go(-p)<5i(p) + - A |A|2 Go{p)Go{-pMp, ~p; 0, 0) 



(21) 



Finally we obtain a closed equation by eliminating F{p) in favor of A{p) through Eq. (jl7p . or more simply through 
the definition Eq. p6)) . Multiplying Ea. (l2T]) by t/(k — p) and summing over p we find: 



6i{k) = - 



^ Vik - p) Go(p)Go(-p)<5i(p) - i A |Ap 



^ Vik - p)Go{p)G„i-pMp, -p- 0, 0) (22) 



If we had only the first term in the right-hand side, we would merely have the linearized BCS gap equation. Accordingly 
we proceed in the standard way [20, 27] to eliminate the bare potential in favor of the scattering amplitude. This is more 
easily understood by proceeding in a formal way. Introducing the operator V , with matrix elements Vk^k' = ^(k — k'), 
and the diagonal operator D with matrix elements Dp^pi = Sp^p' Go{p)Go{~'p) , the T-matrix for two fermions (in the 
center of mass frame) in vacuum satisfies T — V + VDT = V + TDV , where the operator product implies naturally 
the summation on the intermediate variables with = — (27r)^'' J dp dv. This can just be written V^^ = T^^ + D. 
If we note F'{p) the second term in the right-hand side of Eq. (PT|l . we can rewrite Eq. l^^ as Si = VDSi + VF', or 
T>-i(5i D5i + F'. Inserting the above value of V""\ this leads to f^^Si = F' , or 5i = TF'. This is explicitely: 



,{k)^Y.T^{k,p,E)F'{p) 



(23) 



where E is the energy at which the T-matrix is calculated. In the present case, since Gq{p) = [iuj — ek + m]""'^, the 
total energy for the two fermions is E = 2fi {pi for each fermion), as it can be seen directly by integrating the first 
term of Eq. (|22|) over v. 

In standard scattering theory the scattering amplitude /k(k') at energy E = jm = k'"^ /ra (the reduced mass is 
m/2) is obtained from the evaluation of the T-matrix "on the shell", i.e. taking T(k, k',i?) with E — k"^ /m. In the 
present case we deal with a potential with a short range of order tq. This implies that k and k' have to vary on a 
typical scale l/rg in order to have a significant variation of T(k, k',i?). Since in the case of ultracold atoms we are 
dealing with much smaller wavevectors, and we may replace T(k, k', E) by T(0, 0, E). Hence the scattering amplitude 
is merely given by /k(k') = — (TO/47r) T(0, 0, E). On the other hand, in this same regime, the scattering amplitude is 
given by its s-wave component f{k) — — [1/a — \/ —jnE]^^ . This gives T(0,0, i?) = (47r/m)[l/a— ^/ —mE]'^ . 
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Now we evaluate Eq. p3)) for k small compared to I/tq. We see that in this range ^i(k) is essentially a constant, 
which we call simply 5i. On the other hand F'{p) is a rapidly decreasing function [IGJ of p, on a typical scale 
1/a. Hence we may replace in Eq. (|23l) T(k, p, i^) by T(0,0, 2/z), which can be expressed in terms of the scattering 
amplitude. This leads explicitely to: 



m 
An 



- - ^/MHl '5i = J A |A|2 J2 Go{p)Goi-pMp, -p; 0, 0) 



Finally we take into account the fact that the sum over p _appearing in the right-hand side of Eq. 
to the dimer-dimer scattering amplitude T4 in vacuum [iq by: 



T4(0, 0; {0, 2m}) = J2 Goi-p)Go{pM~p,p; 0, 0) 



(24) 

is directly related 
(25) 



whereas: 



Stt 



T4{0,0;{0,-Eb}) = 



27r(2aAf ) 



(26) 



We note that, in Eq. (j25p . we end up with T4 being evaluated at a dimer energy 2/i, while the dimer-dimer scattering 
in Eq. (P5|) is related to T4 evaluated at the binding energy —E}, = —1/ma?. However at this stage the small difference 
between 2/i and — 1/ma^ is unimportant in T^, since it appears in the corrective term which is on the right-hand side 
of Eq. (IM)l . In this term we may also, within the same level of accuracy, use 5i ~ A. This leads finally to: 



■CLM |A| 



(27) 



When we substitute in this expression the relation |Ap ~ STrn/m^a found in section Hill A. we obtain at the same 
level of accuracy: 



1 



M = 



Tin aM 



m 



(28) 



which is the standard mean field correction to the chemical potential, with the proper dimer-dimer scattering length 



IV. COLLECTIVE MODE 



Since the LHY terms come physically from the low energy collective mode, we need first to derive the expression for 
the collective mode propagator, as well as its anomalous counterpart. These propagators can be seen as the extension 
to the superfluid state of the dimer propagator in the normal state. For the normal component of this propagator, 
this is just the full vertex with two entering and two outgoing fermion lines, the two fcrmions having opposite spins. 
In the normal state, in the dilute limit, this vertex is just the single dimer propagator T2 already used in Ref . [l6l| . 
It depends only on the total momentum P and the total energy Q of the two fermions. In contrast with the above 
reference we work in the grand canonical ensemble and the e nerg ies are taken with respect to the chemical potential 
II of single fermions. Hence for all the energies e used in Ref. [16l| . we have to make the substitution e ^ uj + fj,. For 
the single dimer propagator we have to make E n + 2/i. However we keep for simplicity the same notation as in 
Ref. 16] for the various quantities after having made this substitution. Also for coherence, and in contrast with the 
preceding sections, we work in the present section with real frequencies, rather than with Matsubara frequencies. This 
is also more natural physically since we are interested in particular in the (real) frequency of the collective mode. 

Making use of the quadrivector notation P = {P, ft} we have thus, for this single dimer propagator in the dilute 
limit: 



j.(p.J][ 1 ^ 7£^+v/PV4m-0 + 2|/i| 

' m3/2 ^Py4m-n + 2\ii\ m^^ Q ~ 2\fi\ ^ PyAm + Et ^ ' 

where iJf, = 1/raa? is the dimer binding energy, and we have made explicit the fact that /i < 0. Note that physically 
this propagator describes the bound state of the two fermions, through its pole at O = 2|/i| -f P^/4m — E't, , as well as 
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FIG. 5: Bethe-Salpeter equations for the full vertex in the particle-particle channel, a) Normal part, b) Anomalous part. 



the continuum of scattering states where the dimer is broken, through the continuum in the spectral function starting 
at 17 = p2/4m-|-2|/x|. 

In the superfluid we want to obtain the correction to this normal state result to lowest order in |Ap. More precisely 
we want to correct what corresponds to the self-energy if we had a single particle propagator. In the present case this 
is the irreducible vertex in the particle-particle channel, which by definition can not be separated in two disconnected 
parts by cutting two fermion propagators. More precisely, since the bare interaction between fermions is already taken 
into account in T2 (P) , this is the part of this irreducible vertex beyond the bare interaction we are interested in. This 
part Tirriqi, q2', g'l, ^2) ^^'^ irreducible vertex depends on the quadrivectors qi and q2 of the entering particles, as 
well as on those q'l and (73 of the outgoing particles. Assuming for the moment quite generally that the full vertex 
T{pi,p2',Pi,P2) depends also on these four quadrivectors, we can write a Bethe-Salpeter equation, which is shown 
diagrammatically in Fig. [5^: 

^{pi,P2;p'i,P2) = T2{pi +P2 =p'i +P2) + X! ^2(pi +P2 = qi + q2)Tirr{qi,q2;q'i,q'2)^{q'i,q2,p'i,P2) 

+ T2ipi+P2^qi + q2)rt,{qi,q2;~q[,-q2)'i^%-q'i.-q2,P'l,P'2) (30) 

where = */ qdfl / {2tt)'^ . We have made explicit the fact that the normal state dimer propagator T2{P) depends 
only on the sum P = pi + P2 = p'l + p'2 of ttic incoming and outgoing fermions quadrivectors. In Eq. pop the first 
sum is completely analogous to what should be written in the normal state. However we had also to write the second 
sum, which describes the anomalous processes where a dimer can be annihilated into the condensate or can be created 
from the condensate. This leads to introduce, in correspondance with the first sum, an anomalous irreducible vertex 
FjYj, where four fermions are entering (annihilated into the condensate) and none is going out. Correspondingly we 
need to introduce an anomalous full vertex with only outgoing fermion lines. Naturally in order to close our set of 
equations we have to write for this vertex the corresponding Bethe-Salpeter equation, which is shown diagrammatically 
in Fig. [5)3, the major difference with the preceding one being that there is no normal state term: 

r°(-pi,-P2;p'i,P2) = X! ^2(-pi-P2 ^ -qi- q2)Tirr{-q'i,-q2-,-qi,-q2)'r''{-q'i,~q2,p'i,p'2) 

+ ^2(-pi -P2 = -qi - g2)rr„(-'?i,-g2;'7'i,g2)r('7'i:92;Pi>P2) (3i) 

where we have introduced the irreducible vertex F?,^ where there are only outgoing fermions, physically created from 
the condensate. 

Now it is explicit in Eq. (|30p that F(pi,p2;p'i,P2) depends only on P — pi p2. By time reversal symmetry it 
depends also only on p'l + p'2 = P, this last equality resulting from momentum-energy conservation. Hence, just as 
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T2(P), it depends only on the total momentum and energy P, and we merely denote this function as r(P). Similarly 
r"(— Pi, — P2;Pi,P2) depends only on —P = —pi — P2 — v'\ + p'2 ^^^^ '^^ call it merely r°(— P). In this way Eg. ([50]) 
becomes simply: 

r(p) = T2(P) + P2(P) ri„(P) r(p) + T2(p) rf„(p) r"(-p) (32) 

and similarly: 

r'^(-P) = T2(-P) Tirri-P) T'^i-P) + T2{-P) f'}„.{-P) T{P) (33) 
where we have introduced: 

ri„(P)= Y,V,rr{qi,P-qi-.q[,P-q[) (34) 

rf„(P) = rf„(qi, P - gi; -g'l, -P + (35) 

f f,,(P) = ^ ff,,(-gi, -P + gi; (zi, P - gi) (36) 

These equations are naturally quite similar to the ones we had for the single fermion Green's function. Their solutions 
are: 

v(p')^ '^2 ^(~-^) ~ ri.r^(-P) ^ 

[r-'(p)~ri..(p)][T-i(-p)-ri„(-p)]-rf„(p)fr,,(-p) ^ > 

and 

[T2-i(P) - ri„(P)] [r2-i(-p) - ri..(-p)] - r-,(p)ff,,(-p) 

As it could be expected physically, these equations are formally identical to those obtained for pure bosons '25*1 . In 
the same spirit we note, for later use, that we could also introduce an anomalous full vertex F" with only ingoing 
fermion lines. It is easily shown, by writing the equivalent of Eg. ([5^ with this vertex, to be related to the preceding 
ones by f '^(P) f f^^(-P) = r'^(-P) rf„(P). This gives: 

(P) 

r'^fP) = "'''^ ' ^ (39) 

[T^\P)-T,rr{P)][T2\~P)-T-„.r{-P)]-Tf„{P)ff„{-P) ^ ' 

We have now to find the expressions of Tirr, T?.^ and F?.^ within our expansion in powers of A. Since A corresponds 
to diagrams with two ingoing fermions and A* to two outgoing fermions, particle conservation implies that the 
expansions start at second order with ri„. proportional to AA*, F?,^ proportional to A^ and f?^ proportional to 
A*^. Since it can be seen that it is not required to go to higher order, the expressions of our irreducible vertices 
will be given by these second order expressions. Since we have written explicitely the various factors A and A*, the 
coefficients in the expressions of Firr,F?,^ and Tf^^ correspond to diagrams containing only normal state propagators. 

For Firr (see Fig. ^^), we have to write the expression for all the normal state diagrams with one ingoing dimer 
propagator and one outgoing dimer (from the very definition of T{rr)- Moreover there are in addition two lines 
ingoing from A* and two lines outgoing to A. These last two couples of lines have to be treated as dimer lines, which 
means that they are not allowed to interact immediately in the normal state diagrams we are drawing, because such 
interactions are already taken into account in A and A*, and allowing them would amount to double counting. In 
summary we have to write all possible normal state diagrams with two ingoing dimer lines and two outgoing dimer 
lines (these diagrams are automatically irreducible in the sense of the definition of Fi^r, that is with respect to a single 
dimer propagator, since the four ingoing normal state lines can not disappear). This is depicted diagrammatically in 
Fig. [5^. But this set of diagrams is just by definition the normal state dimer-dimer scattering vertex T-^ introduced in 
[l6(|. More precisely, since the total momentum-energy linked to A and A* is zero and the total momentum-energy of 
the ingoing lines is P, this vertex is T4(P/2, P/2; P/2) with the notations of Ref.il6] which gives for Fj^r the explicit 
expression: 

T,„{P)^\1^?tJ^,^-^\ (40) 
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FIG. 6: Diagrammatic representation for the irreducible vertices coming in the collective mode propagator, at second order in 
A. a) Normal irreducible vertex FirriP)- b) and c) Anomalous irreducible vertices rj?,r(P) and f?.^(P). 



With completely analogous arguments for F?^, we obtain: 



(41) 



and similarly: 



(42) 



corresponding to the diagrams found in Fig. and Fig. \Et- In Eq. ((4T|l and (|42ll we had to put a 1/2 topological 
factor corresponding to the fact that, for example, exchanging the two factors A in Eq. (|4T|l (that is precisely the set 
of diagrams they correspond to) does not produce a different diagram, whereas this diagram would appear twice if 
we had just written A^ T4. 

Finally it is worth pointing out that T4 satisfies T4{—pi,p2',P) — Tji{pi,p2', P)- This is physically obvious since it 
corresponds to state that exchanging the arguments of the two ingoing dimer lines does not change the value of the 
vertex, since we exchange two bosons. But this property can also be proved directly from the integral equation [l6j 
satisfied by Tj^. Similarly we have T^ipi, —p2] P) = Ta{pitP2] P)- 

Let us consider now more specifically the collective mode dispersion relation. It is given by the pole of the collective 
mode propagator, that is: 



[T-\P) - T.rriP)] [T^'i-P) - T.rri-P)] " Tf,, (P)f f,,, (-P) = 



(43) 



Naturally we have first to check that for zero wave vector P = 0, the frequency Vl is zero, that is P = is solution of 
this dispersion relation. For this case we have from Eq. (|29p for the bare dimer propagator: 



7.3/2 



47r 



On the other hand we have 



ri„(0) 2F«,(0) 2f«,(0) 



T4(0,0;0) = — mVflM 
Idtt 



(44) 



(45) 



where the last equality comes from Ref . . and we have used the fact that, in r4 (0, 0; {0, — £^6 + 2|/i|}), we may, at 
the order we are working, replace by its zeroth order value, namely |/Li| = Eb/2 (we recall that we have to shift all 
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the frequencies of Ref. [ig'I by 2/i). Since we have already found that •\/2|/i| = \^Eb — |Apm'^/^a^ aM/8, we can easily 
check that Eg. (1^5]) is indeed satisfied in this case. 

Now, as we will see, we are more specifically interested in the collective mode dispersion relation for wave vectors 
|P| ^ l/a and frequencies Q <ti Eh. Since and 2|/i| are of the same order of magnitude, we can expand the bare 
dimer propagator to first order as: 



T,^\P) = T,-\0) —= —-n] (46) 




On the other hand the typical scale of variation of T4 on its arguments is l/a for the wave vectors and Eb for the 
energies, since we deal with a normal state quantity and there are accordingly no other scales in the corresponding 
problem. Hence, for the wave vectors and frequencies we are considering, we can still take: 

On the other hand the above result at P = implies also that we have T^^iO) - ri„(0) = -(l/2)rirr(0) = -|r° ^(0)|, 
and Eq. for the collective mode frequency flcmi'P) can be rewritten as: 

\ 4m / Am \ 4to / to^ 

where we have set fiB — 'i''^\/^\jj]^irr{0)/'fn^^'^ and in the last step we have used the lowest order expressions for 
and |Ap to calculate this quantity. Naturally Eq. (|^5)) is just the Bogoliubov dispersion relation, with the proper 
scattering dimer-dimer scattering length. 



V. COLLECTIVE MODE CONTRIBUTIONS TO THE SELF-ENERGIES 



As we have already indicated in the introduction, the LHY term arises from the collective mode contributions to 
the self-energies. If we were to continue naturally our above procedure and go to next order in A, i.e. to order A^ for 
example in the calculation of the normal self-energy, we would find singular results (as it can be checked by proceeding 
to such a formal expansion in the explicit calculations performed below). This is due to the fact that the collective 
mode is gapless, and the lack of gap in the excitation spectrum leads to a singularity in the perturbative expansion. 
Hence, instead of the next order terms, we include the whole contribution coming from the collective mode. Avoiding 
in this way to proceed immediately to a A expansion is equivalent to a series resummation, which allows to take care 
of singularities. 

Writing the collective mode contributions to the self-energy is formally easy, when one notices for example the 
analogy of the normal part of the collective mode propagator with the product of A* and A. Indeed the term A* acts 
as a 'source' of fermions in our diagrams, coming physically out of the condensate, while A acts as a 'sink' of fermions 
going into the condensate. These source and sink related to the condensate are required since, at T = 0, < 0, 
no fermions are present except coming from the superfluid. However when we allow to insert the normal part of a 
collective mode propagator, the starting part of the propagator acts as a sink for the rest of the diagram while the 
end part acts as a source. Hence we have just to substitute AA* by the normal part of a collective mode propagator. 
We will proceed in a similar way for the anomalous parts. Naturally we will have to take care that the self-energy 
contributions already calculated above are not counted twice, by removing the corresponding terms. 



A. Contribution to the normal self-energy 



Hence the collective mode contribution 'Ecmip) to the normal self-energy is obtained by replacing, in our lowest 
order expression, the product AA* by the collective mode propagator, more specifically its diagonal component r(P) 
given by Eq. (P7|) . as depicted diagrammatically in Fig[71 We have naturally to sum over the variable P of the collective 
mode. Corresponding to Eq.®, this leads to: 



Sc™(p) = ^r3(p,p;p + P)r(P) 
p 



(49) 
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FIG. 7: Diagrammatic representation of Ecm(p)- The diagonal component r(P) of the collective mode propagator is depicted 
by the heavy line 
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FIG. 8: First order term in |Ap in the collective mode contribution to the normal self-energy Ea. (|49|l depicted in Fig[7] 



In principle this expression is improper because it contains, when we perform its expansion in powers of |Ap, 
contributions of zeroth and first order in |Ap which have already been taken explicitely into account in our above 
expansion. However in the present case these two contributions are zero. For the zeroth order term, this is fairly 
obvious physically since it corresponds to a normal state diagram with a normal state dimer propagator, and at T = 
no dimers are present which makes all diagrams of this kind equal to zero. Mathematically this would happen because, 
in this diagram, the presence of a normal state dimer would necessarily imply, since there is no down-spin entering 
line, the presence of a down-spin loop (if we calculate the up-spin self-energy) which can not be present as explained 
in Appendix [Bl 

The conclusion is the same for the first order term, but the argument is slightly more complicated. The factor of | Ap 
in this term is J^p T3{p,p;p+P)[T2{P)]'^T4{P/2, P/2; P/2), which is a normal state diagram (see FiglS]). It is actually a 
part of the diagrams belonging to T3(p,p;p) since it has one entering fermion line and one entering dimer line, with the 
same for the outgoing lines. However we have seen, as it is detailed in Appendix[Ul that only the Born term in Tt, gives a 
non zero contribution to the particle number. On the other hand the term J^p Tsip^ P'-, P+ P)[T2(P)]'^T4{P / 2, P/2; P/2) 
does obviously not contain the Born term of T3 since this last one reduces to a simple propagator, whereas the presence 
of the factor T2(P) implies at least two lines, with interaction between them. Actually one can see more precisely 
that our term is equal to zero, since it contains necessarily propagators going backward in time (see Appendix [Bl and 
p. Taking the Born contribution for T3 and T4 in J2p T3{p,p;p + P)[T2iP)]^Ti{P/2, P/2; P/2) gives a good example 
of this. Finally, after these considerations, we see that we can still use Eq. (P5)) for Ecm(fc)- The conclusion will be 
different for the mode contribution to the anomalous self-energy discussed below. 

Now when we look for the contribution ricm to the density corresponding to Eq. (j49p . we have to replace G{p) in 
Eg. pU)) by Go(p)I]cm(p)G(p), just as it is done below Ea. (fTU|) . Quite similarly it will be enough to replace here G{p) 
by Go{p) to obtain the lowest order contribution of the collective mode we are interested in. This gives: 



= X! ^3{p,p;p 
p,p 



P)T{P) [Go{p)f 



(50) 



where again J2p = ~(27r) ^ / dpdv and p = (p, iv). Now the argument runs just as below Eq. ([TT)) . We integrate first 
on the frequency variable v of the quadrivector p. As shown in Appendix [Cl it turns out that, except for the Born 
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contribution, T^{p,p;p+P) is analytical in the upper complex half-plane of the variable v. Since [G'o(p)] = [ji' — Cp]~^ 
is also analytical in this half-plane, one obtains by closing the integration contour at infinity in this half-plane that 
only the Born part —Go{P — p) of T3(p, p;p + P) contributes to the result. Hence for the integration over i/, we have 
just to calculate: 



dv 



1 



1 



1 



and we are left with: 



in 



dpdP 



[iQ - 2\fi\ - Ep - ep-p]^ 



dn 



r(P) 



(27r)6 27r [tfl - 2\fi\ - - ep_p]2 



(51) 



(52) 



The factor of T{P) is analytical in the upper complex half-plane of fl (which corresponds to the negative frequency 
half-plane when we go back to the standard frequency language). Hence by closing the fl contour in this half-plane, 
we are left to consider the singularities of T{P). Physically they correspond to the two-fermions excitations of the 
system, that is the collective mode and the broken molecule excited states. In the standard frequency language, the 
corresponding singularities are on the real frequency axis (the imaginary axis for Q), on the negative as well as on the 
positive side. 

Let us first consider the low frequency singularities. For a fixed value of P, we have just as singularity a single 
pole corresponding to the collective mode frequency. Let us call V,cm{P) the positive solution of Eq. (|48|) . Our pole is 
located at — if2cm(P)- In Eq. ((37| . which gives r(P), for low frequency and wavevector, the denominator becomes 
m^{nl^(P) - r22)/(i287r2|/i|), while from Eq.® the numerator is -[(l/2)ri^^(0) + m^/'^(P^/im + n)/{8TTy/2]iI\)]. 
After going to imaginary frequency f2 —^ zf2, this gives in T{P), for the pole at iflcm(P)-, a residue: 



Z{P) 



I m 



3/2 



1 



4m 



MB 



^cm (P) 



On the other hand, for small VL and P, the p integration in Eq. ([52|) reduces to: 

dp 1 m^/^ 



This leads to: 




(2^)3 [2|Ai| + 2ep]2 87rV2R 



zZ{P) 



47r2 



dx X 



+ 1 



(53) 



(54) 



(55) 



where in the last step we have made the change of variables P 
our final result is: 



{Ara^BY^'^x. The numerical integral is a/2/3 and 



37r2 



37r2 



(56) 



where in the last equality we have made use of the lowest order result for |A|2, given by Eq. (fT3|) . Naturally, once we 
have taken care of Ta, we are essentially back to the elementary boson theory 25, 28] and one recognizes in Eq. ([551) 
the standard "depletion of the condensate" for the dilute Bose gas, with ^ = (4m/^s)~^^^ being its coherence length. 

We have first considered the contribution of the low frequency singularities. We have now to finish up by showing 
that the higher frequency ones give a negligible contribution. First we see that, in Ea. (j52p . the denominator plays in 
this respect an unimportant role, since its modulus is larger than (2|/^| -|- ep)^. Hence to have an upper bound for the 
contributions we are investigating, we can replace this denominator by this lower bound, which gives after integration 
over p an unimportant factor irt'^'^ / {A'it\^\^/'^). Hence we are again left with the evaluation of J d^T{P). 

It is now more convenient to go back to real frequencies (equivalent to the change of variables = to). The 
integration is now over the imaginary uo axis and is closed in the negative lo half-plane. The contour can be deformed 
to merely enclose the negative oj axis, and corresponding to the singularities on this axis, Imr(P) has a jump across 
this axis. Hence we have: 



-OO 7 

— lmT{{V,uj + ie]) = 



duj lmT~^ {{P, LJ + 
~ \T-^{{P,oj + ieW 



(57) 
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FIG. 9: Diagrammatic representation of Acm(p)- The diagonal component r(P) of the collective mode propagator is depicted 
by the heavy line with a single arrow and the off-diagonal one r"(P) by the heavy line with two arrows. For the two diagrams 
the slash in the heavy lines means that the terms which are first order in |A|^ have been removed. 

Now, from Eq. ([57|) . r~^({P, u; + ie}) contains a first term T^^{P) — Tirr{P) which has aU its singularities on the real 
positive uj axis, corresponding to the continuous spectrum of the broken dimer (for T2(P)) and of the broken pair 
of dimers (for TirriP) = jAp (P/2, P/2; P/2)). Hence it is real on the negative CO axis and does not contribute 
to ImF"! in Eq.dlTl). From Eq.dSZ]) the second term of is -Tf^^{P)rf^^{-P)/{T^^{-P) - ri„(-P)), which is 
directly proportional to |A|* from the product r?,^(P)r?,^(— P). Accordingly we have proved basically our point since 
our result, being proportional to lAj** ~ n^, is indeed negligible compared to the result Ea. (j56p . which is proportional 

to 71^/2. 

Actually, to complete our proof, we have to be more careful and evaluate the coefficient of | A|^. This is quite obvious 
since Eq. (j56p is in fact a contribution contained in Eq. (j57p (we have made no approximation to reach Eq. (j57p ). This 
detailed analysis is performed in Appendix ID] 



B. Contribution to tlie anomalous self-energy 

We will now proceed in a similar way to obtain the mode contribution Acm{k) to the anomalous self-energy A(fc). 
This contribution is obtained by replacing, in the diagrams which lead to the second term of Eq. (|17p . any second order 
product AA*, A A or A* A*, by the corresponding collective mode propagator. Accordingly we replace specifically 
AA* by the diagonal component r(P), and AA by the off-diagonal one f°(P). This is depicted diagrammatically 
in Figini Just as in Eq. p?|) this last term comes with a topological factor 1/2 to avoid double counting. On the 
other hand, since the two dimer outgoing lines in the first term play inequivalent roles, such a topological factor does 
not apply to this first term. It is worth noticing that, since the collective mode propagator depends only on the 
total momentum-energy four-vector, we do not meet at this stage the problems linked to the fc-dependence of A(fc), 
discussed in Appendix [X] 

Just as in section ITlI Bl we should not include the terms of $(gi, g2;P2, P) leading to reducible contributions to the 
self-energy. These are easily obtained from the diagrams displayed in Fig[3]by making the above substitution of AA* 
by r(P). This implies in the same way the replacement of ^{qi, q2;p2, P) by ^' {qi, q2;p2, P)- However, as in section 
IIIIBi it is more convenient to add these reducible terms to both members of the equation for A. In the left-hand 
side one obtains [Go{p)]^^F{p)[Go{—p)]^^, while in the right-hand side one recovers instead of Finally we 
restrict ourselves to the contribution of the low frequency collective mode, the higher frequency contributions being 
expected to be of order A'"^ and thus negligible, just as in the preceding section IV Al Accordingly the variable P in 
r(P) and r°(P) is small. This implies that $(p, —p; P/2, P/2), which comes as a factor of r(P), can be replaced by 
^{p, —p;0,0) at the order we are working. Similarly $(p, — p;P, 0), which comes as a factor of T°'{P), can also be 
replaced by $(p, — p;0, 0). After multiphcation by Go{p)Go{~p) this leads to the equation: 



Fip) = Go{p)Goi~p)S,{p) + AGo{p)Goi-p)Hp,-p;0,0) 



1 



A\' + J2riP) + lY.r{P) 



(58) 



which comes in place of Eq . (|T9|) . 

In the preceding section IV Al we have already calculated X]pr(P) at the level of Eq. ([55|l and found J2p^{P) = 
8ny/2\fi\ rn~^/2 Ucm: where ricm is defined by Eg. ([56)1 . In principle the calculation of f °(P) follows the same line 
as the one of J^p r(P)- One performs first the frequency integration by noticing that, for fixed P, r°(P) has a pole 
a,t — iilcmCP)- The corresponding residue Z°(P) is obtained from the residue Z{P) for r(P) by comparing Eq. p8|) 
with Eq.jifl). One has to multiply Z(P) by F? ^(P)/ [T2"^(-P) - Fi^^(-P)] , evaluated at the pole, to obtain Z'^(P). 
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This gives merely: 



(59) 



Now if we proceed as in Eq. (|55p we come to an integral with an ultraviolet divergence, which would invalidate our 
procedure. 

The reason for this problem is that wc have not yet consistently removed in our calculation the terms of order zero 
and one in |Ap, as we have done precedingly. That we should do it was already appearent in Eq. (|58p by the notation 
f'" instead of (which implies also a multiplication by A* /A). We recall that these terms have to be subtracted 
out because they have already been taken into account in our preceding lowest order calculation. However in the 
above preceding cases, it turned out that these terms were equal to zero, which made the subtraction transparent. 
Here this is not the case. The most efficient way to perform this subtraction is to note that, from Eq. ([39|) . r'^(P) 
contains already an explicit factor r?,^(P) in the numerator, proportional to |Ap, which gives directly the factor [Ib 
in Eg. ((59)) . Hence the term of order |Ap to be subtracted is merely obtained by taking the other factors to order 
zero in |Ap, that is by replacing 17cm (P) in Eq. ([59|l by its zeroth order expression P^/4m. Actually this argument is 
a bit too fast. It misses the fact that we may at the same time evaluate $(p, — p; 0, 0) in Eg. ([55]) with the chemical 
potential taken as —Eb/2, even for the (l/2)|Ap term, as we had done in Eg. ((25)l . Hence this whole step is taken up 
in more details in Appendix [Ej We find in this way: 




1 Am 



(60) 



Making again the change of variables P = (4m/iB)^/^x, this gives 

dx , ^ — 1 



\/a;2 + 2 



TT 



7-3/2 



(61) 



Hence the last bracket in Eq.dHS]) can be written (l/2)|Ap + (5/2) J2p ^(P) 



The rest of the argument proceeds just as in section UlIBl provided (l/2)|Ap is replaced by the above bracket. 
This leads to the replacement of Eq. (^7)) by 



- - ^/2m\fi\ 



207rV2R 

™3/2 



(62) 



We may eliminate |A| by making use of the number equation, which reads now from Eq.([T 



This gives: 



m^|Ap 
87r(2m|M|)i/2 



(63) 



(64) 



where in the last equality our replacement of 2m|/i| by 1/a^ introduces only errors of order in the equation of state. 
To the same order this leads us to: 



2/i = -- 



1 



2T:aM 



n + 



(65) 



In the last term we may replace ricm by its expression Eq. (|56p ricm = (47rriaM)'^^^/37r^. This gives the familiar 
Lee-Huang- Yang result for the chemical potential of a dimer, taking its binding energy as reference: 



2n + 



1 



Anna 
2m 



M 



32 / na 



3 \ 1/2- 



(66) 
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VI. CONCLUSION 

In this paper we have presented a general theoretical framework to deal with fermionic superfluidity on the BEC 
side of the BEC-BCS crossover. Our approach involves no approximation in principle, but in order to obtain explicit 
results we have to proceed to an expansion in powers ot the anomalous self-energy, which is basicalUy equivalent to 
an expansion in powers of the gas parameter na^. We have applied our method to the T = thermodynamics and 
calculated the first terms in the expansion of the chemical potential in powers of the density. We have in this way 
shown explicitely that the mean-field contribution is indeed given by the standard expression in terms of the dimer- 
dimer scattering length gm- We have then proved that the Lee-Huang- Yang contribution to the chemical potential 
retains also its standard expression from elementary boson superfluid theory provided om is used in its expression. 
Departures from elementary boson theory appear at the order n^. Naturally our approach can be used to obtain other 
physical quantities and extended for example to finite temperature and dynamical properties. This will be considered 
elsewhere. 
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APPENDIX A: EXPRESSION OF E'^' IN TERMS OF Tg IN THE GENERAL CASE 

If we note T^°'{p,p' ,p") the sum of all normal diagrams, with one single atom p line going in and out, two entering 
atom lines (p',T) ^nd (— p', i) corresponding to an ingoing dimer, and two outgoing atom lines (p",T) s-^d {~p" 1 1) 
corresponding to an outgoing dimer, we have: 

= j2 ^*{p'W)Tr{p.p\p") (Ai) 
p' p" 

In order to obtain an expression for T^°'{p,p' ,p"), one can find in the spirit of Ref. 16] integral equations for this 
quantity, and then identify in them a quantity which is precisely T3(pi,p2; P). Another slightly more direct way is to 
separate, in Tg °, all the diagrams where or 1 dimer lines are present, right at the beginning or right at the end of 
T^°'{p,p' ,p"). We do not enter into the details but only give the result, on which the above mentionned separation is 
clear: 

T^%p,p',p") = -G{-p) Sp^p, Sp,p„ + [G{~p')f G{p') T2{p ~ p') V.p" 

-HG(p')G(-p')G(p")G(-/)r2(p-p')T2(p-p")r3(p',p";p) (A2) 

where Spy is the Kronecker symbol, and T^Op) is the normal state dimer propagator. One can check, by making use 
of Eq. (jA2|) and of the integral equation [lg| satisfied by T3, that Eq. (IAl[) reduces to Eq. ^ when A{p) does not 
depend on p. 

APPENDIX B: ABSENCE OF LOOPS IN NORMAL STATE DIAGRAMS 

Since we are at T = with a negative chemical potential, there is no Fermi sea and no possibility of fermionic holes. 
When we consider the propagation of a fermionic atom, it has to be created before being destroyed. This corresponds 
to the fact that the free fermion propagator is retarded. In (imaginary) time and wavevector representation we have 
for the free atom Green's function: 

Go(p,t) = -e--"0(T) (Bl) 

where 9(t) is the Heaviside function. 

We consider now the contribution coming from any loop in a diagram. For clarity we do not write the wave vectors, 
which do not play any role in the argument. The loop starts for example at a vertex with time tq, with successive 
interactions at times ri, T2, r„. This implies in the term corresponding to this diagram the presence of the product 
Go(ro — ti)Go{ti — r2)....Go(Tn — tq). In order to have each term being non zero in this product of retarded Green's 
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function, we need tq > ti > T2 > ... > t„ > tq, which is impossible to satisfy. Hence one of the Green's function in 
this product will necessarily be zero, and the contribution of any diagram with a loop will also be zero. Note that 
the argument applies also if there is a single Green's function in the loop If we were to carry out calculations in 
frequency, rather than time, representation the corresponding free atom Green's function is analytical in the upper 
complex half-plane, and the above result would appear through integration over frequencies with contours closed in 
the upper complex half-plane. 

We have considered for consistency the imaginary time formalism. Naturally exactly the same arguments apply if 
we consider the real time representation. 



APPENDIX C: ANALYTICAL PROPERTIES OF Ts{k,k;k + P) 

We give first an explicit proof of the analytical property of T3 indicated in the text. We start from the integral 
equation satisfied by T3 with at first essentially the notations of Ref. : 

T3ip,,p2;P) = -Go(P - pi - P2) - 5] Go(P - pi - q)Go{q) T^iP - q) niq,p2; P) (CI) 

q 

We follow Ref. 16] and perform the frequency integration in the second term, making use of the fact that the only 
singularity in the lower half complex plane comes from Go{q) and is located at q^/2m — iO_|_. We write then explicitely 
the resulting equation for pi — K — {K,in}, p2 = k = {k, iw}, while P = {F,iv} is replaced by fc -I- P, keeping in 
mind that the frequencies pi and p2 have to be shifted by (and fc + P by 3/i) and that we work now with imaginary 
frequencies. This gives: 

1 47r /■ dq 1 



1 47r /■ 

r3(i^, fc; fc + P) - - ^ ^ z(i/ ~n)-{P- K)2/2m " ^ J ' 



fi + i{i'~n)-{P - K)2/2m m3/2 J {2tt)^ 2fi + ~ n + uj) -q2/2m-(k + P - K - q)2/2TO 

ni{q,qy2m-f,},k;k + P) (C2) 



4- ^(P + fc - g)2/4m - + oj) + 3/i - q2/2TO) ^ 



i{v + uj) + 3/i - q2/2m - (P + fc - q)'^/Am + Et, 

When we write this equation for K = k, we see that the first term, which is the Born contribution, has a pole in the 
upper complex half-plane, at lu — v + i(|/i| -I- (P — k)2/2m). On the other hand, in the integral, the first factor does 
not depend on frequency uo. The singularities coming from the dimer term will be in the lower complex half-plane 
withlmw < — (3|/i| — -Eh) < (since we will have ~i?f,/2). Regarding the r3({q, q2/2m — /i}, fc; fc -I- P) term, if we 
write Eq. (jC2p for the variable K evaluated on the shell, i.e. K = {K,K2/2m — /x}, we obtain an integral equation 
for Tgdq, q2/2rn — /x}, fc; fc -(- P). In this equation the Born term does not depend on frequency. The first term in 
the integral has its poles satisfying Imo; < — 3|/i| < and the dimer term has again its singularities in the lower 
complex half- plane. Hence by a recursive argument (or equivalently by iteration of the equation), we conclude that 
the singularities of P3({q, q^/2m — /i}, fc; fc + P) are also in the lower complex half-plane. This shows that indeed, 
except for the Born contribution, r3(fc, k;k + P) is analytical in the upper complex half-plane of the frequency variable 

UJ. 

Naturally we do not expect such a result to come by chance and we discuss now its physical origin. It is more easily 
seen in real time, rather than frequency, representation. The particle number Eq. ljlOp is merely given by: 

E( 

k 



n= > (cI^(e)cfeT(O)) (C3) 



with e ^ 0+. The average corresponds to a full propagator Gfcf(— e) in time representation, where the particle is 
created at time t = e, later than the time t = at which it is destroyed. In standard [2^ diagrammatic representation 
(see Fig llOl a)), with time going from left to right and the propagator starting at the destruction time and ending at 
the creation time [25j , we have to draw diagrams where the "f particle propagator starts from t = and goes globally 
forward (that is to the right in Fig llOp in time, ending up at t = e. However when we write the perturbation expansion 
of this propagator, we have to use free particle propagators Go which are retarded and go only forward in time (i.e. 
the arrows go to the left in Fig [TUl when we use the standard Green's function representation [25]). Hence the zcroth 
order in our perturbation expansion is obviously zero. When we go to second order, corresponding to Eq. ljlip . we 
can use A* which creates, at some given time, a pair of particles (q tj~Q i) (this process is instantaneous in time 
because, at lowest order, A does not depend on frequency). Similarly A annihilates a pair of particles. The simplest 
combination of these ingredients, giving a non zero result for n, is shown in Fig llOl b): A* is located at positive time, 
and A at negative time, allowing to have Go's going only to the left. This is just the Born term of r3. However one 
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FIG. 10: a) Diagrammatic representation of the propagator Gfe|( — e) involved in the calculation of the particle number n. 
We use the standard representation where the propagator starts at the destruction time and ends at the creation time, b) 
Born term in the diagrammatic expansion of the propagator represented in a), c) Example of a diagram which gives a zero 
contribution due to the presence of a loop. 



sees easily that it is not possible to draw any other second order diagrams having the same property shown in Fig llOl 
b). They would necessarily imply Go's going at some stage backward in time (i.e. with arrows pointing to the right 
in Fig fTU]) . which is not allowed. Indeed we have to dispose of the t propagator starting at t = 0, and this can only 
be done by annihilating it with A at some negative time. Similarly the ] propagator ending up at i = e has to come 
from A*, acting at some positive time. Then we are left with a | propagator which can only go from A* to A. Any 
interaction between the t and J, propagators is forbidden since this would merely amount to double count processes 
which are already taken into account in A* and A. In this way wc end up with the conclusion that only the Born 
term leads to a non zero contribution to n. 



APPENDIX D: DETAILED ANALYSIS OF THE NEXT ORDER CONTRIBUTION TO THE NORMAL 

SELF-ENERGY 

In calculating the coefficient of |A|'' in Ea. ((57)) . it is clear that we can evaluate all the other factors to zeroth order 
in A, since we have already an explicit factor |A|*. This implies in particular to make 2|/i| = _Eb = \/{ma?) in these 
terms. This leaves us with: 

Trnp) = — [Tt„{P)f'}„{-p)T^{-p)] |r2(p)p (Di) 

with again P = {PjCj + ie\. Naturally this result could have been obtained more directly, by a straight expansion 
of r(P) in powers of A, instead of resumming a whole perturbation series, as we have done in Section [IVI and then 
expanding the result. 

We may first apply this formula for the collective mode contribution. In the present case it comes from the 
factor T2{—P), which has a pole for uj = — P^/4to. This is indeed the frequency of the Bogoliubov mode, 
once the frequency is high enough for the zeroth order approximation in A to be valid. Correspondingly we 
have lmT2{—P) = Stt'^ / {m'^ a) 6 (oj + P^/4to). On the other hand, we have for this value of the frequency 
T2{P) = (47ra/m)[l — y'l + 'P'^a^ /2]~^ . This has to be inserted in Ea. (|Dip and the integration over P has still 
to be performed. However if we concentrate on the low P range (where Tf^^{P)Tf^^{P) ~ jAI** (0, 0; 0) /4), we 



22 



obtain a divergent integral / dP/P^ / d|P|/|Pp. This is due to the fact that, in this range, our zeroth order 
handling in |A| is not careful enough. From Eq. pS]) . we see that the dispersion relation has to be modified for 
|P|m (w/^s)^/^ ~ (najv/)^^^ |A|. If we evaluate the low frequency contribution by putting a lower cut-off |P|m in 
the integral, we obtain J d\P\/\P\'^ ~ 1/|P|,„ ~ n~^^'^ ~ |A|~^/^ and the overall result is of order n^^^. This is just 
the contribution Eq. (|56p . where we have obtained the precise coefficient by our careful handling of this low frequency 
domain. Now, if we define the low momentum region as |P| < Pq, where Pq is large compared to |P|m, but small 
compared to 1/a, our above treatment of the low frequency region is fully valid in this domain. On the other hand 
the integral in Ea. (j55p has already converged at the upper bound xq — Po/(4™/iB)^^^, and we may replace it by oo, 
just as we have done in Ea. ((55|) . Hence the contribution from the domain |P| < Pq is exactly given by Eg. ([56]) . On 
the other hand, there are no divergences in the domain |P| > Pq and its contribution is of order | A|^, and accordingly 
neghgible compared to Eg. ([55]) . 

To be complete we have also to verify that there are no other sources of divergence in Eq. (|Dl[) . analogous to what 
we have found at low momentum. Since there is no singularity arising at finite frequency or momentum, the only 
possibility comes from the large frequency \u!\ > f2c, large momentum |P| > Pc region. In order to evaluate the 
contribution in this domain, it is more convenient to go back to the imaginary axis contour by performing backward 
the first step of Eg. (|57)l and write this overall contribution as: 

'/ '^Tf„{p)nA~p)n{-p)T,{p) (D2) 

The scattering amplitudes T4 (P, 0;0) and T4 (0, — P;0), introduced by r?^(P) and f?^(— P), may be shown to be 
equal and are dominated as expected in this range by the Born contribution. They are real on this imaginary 
frequency axis. The Born term has been obtained explicitely in [l^, and T4 (P, 0;0) is proportional to l/(|P||r2|) for 

large |P| and |r2|. From Eg. ((29l) there is in the modulus of the integrand a factor 1/ [fP + (P^/4m)^] ^^'^ coming from 
r|(— P)T2(P)- Then the overall integral is easily seen to be convergent for large |P| and \ which completes our 
proof. 



APPENDIX E: DETAILS ON THE CONTRIBUTION OF THE COLLECTIVE MODE TO THE 

ANOMALOUS SELF-ENERGY 

Since there is no zeroth order term in powers of |A| in r°(P), we have only to remove the term proportional to |Ap 
to avoid double counting, and go in this way from r'^(P) to f (P). The corresponding diagram in F{p) is shown in 
FigHH This leads to: 

r(P) - r^iP) = -Tt„{P)T2{P)T2{~P) (El) 

where r?,^(P) is given by Eq. (|¥T|) . Let us first replace the dimer propagator T2{P), given by Eg. ([^5]) . by its value 
r2 (P) where /i is replaced by its zeroth order value ~Ei,/2. When we consider the contribution from the low frequency 
and momentum contribution, we may replace r?,^(P) by r?,^(0), and the pole of (— P) at = iP^/(4m) gives: 

- JZ f l-K.JP)n'(P)TS(-P)] - (E2) 

This is indeed the second term in the right-hand side of Eg. ([SOI) since in this last equation, which turns out to be of 
order |Ap, we have to make consistently 2\fi\ — Pf, (and similarly we have to take /is = 47r-\/Pbrirr(0)/m'^/^, with 

ri„(o) - 2(AVA)r«,,(o)). 

We are now left with the evaluation of the difference between our original expression and the one we have just 
considered, namely: 

- rt„iP) [r2(P)T2(-P) - T^{P)T^{^P)] (E3) 
which gives at the level of Eq. , after insertion at the level of Eq. ([55]) and summation over p and P: 

- Ia* A2 J2 Ti (0, P; 0) [r2(P)T2(-P) - Tl{P)T^{-P)\ Ti{P, 0; 0) (E4) 
p 

where we have used Eg. (|4T|) for r?^(P) and the relation (see Ref. [lB|) between r4 and namely r4(0, P;0) = 
V Go(-p)Go(p)<i>(-p,p;P,0). 
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FIG. 11: Diagrammatic representation of the difference between ^"(P) and r''(P) 

We show now that this difference is just what is required to be allowed to replace ^ A | ApT4(0, 0; {0, 2/i}) = 
i A I Ap Y.p Go{-p)Go{p)^-p,p; 0, 0) in Eq.JSSl) by i A | A|2T4(0, 0; {0, -Eb}), which is directly related to the dimer- 
dimer scattering scattering length qm by Eg. ((26)) . Indeed it has been shown in [l^ that scattering amplitude 
T4{pi,p2', P) satisfies quite generally an integral equation analogous to the standard integral equation for elemen- 
tary bosons, namely: 

T4(pi,p2; P) = Ai{p,,p2;P) + ^Yl ^4(pi, g; P)T2{P + q)T2{P - q)T^{q,p2;P) (E5) 

where A4{pi,p2', P) is the sum of all irreducible diagrams with respect to cutting two dimer propagators, which is 
itself obtained from an integral equation with explicit kernel p^ . Formally this equation can be rewritten in terms 
of operators as: 

AJ 1 = T^^ + D (E6) 

where we have introduced the diagonal operator D{q; P) = {l/2)T2{P+q)T2{P—q) (the factor 1/2 is again topological). 
We consider now this equation for pi ~ P2 — and P = 0. With respect to the energy variable corresponding to P, 
we take the two values 2/i and —Eh (we note the operators with this last value with a superscript 0). The crucial 
point is that the difference between A4 and A4 is of order 2/^ + Ei,, that is the expansion of A4 in powers of 2/i + Ei, 
is regular. Such a difference is negligible at the order we are working at. The regularity of the expansion can be seen 
from the integral equation (Eq.(B2) of Ref . [l6j ) satisfied by A4, which is completely analogous to the one satisfied by 
$ (Eq.(9) of Ref.[16|). except that the last term (with two dimer propagators) in the equation for $ is not present for 
A4. This integral equation can be transformed in the same way, leading to more explicit expressions (see Eq.(ll) and 
Eq.(13) of Ref.flBl). On these expressions one checks easily that, for fi in the vicinity of —Ei,/2, the propagator T2 is 
evaluated for values of its argument which are never in the vicinity of its pole. Hence one can proceed to a regular 
expansion. By contrast one sees also on the equations for $ that, for the term containing two dimer propagators the 
argument does reach the pole value, which makes the expansion singular. This explains why we can not argue that 
the difference between T4(0, 0; {0, 2/i}) and r4(0, 0; {0, — _Ef,}) is of order 2/i + Ei, and negligible. The equation for 
P = {P, -Eb} reads: 

[Al]-' ^ [T^]-' + D'' (E7) 
Writing the difference with Eq. (|E6|) and taking into account A4 ~ A4, we have: 

[T°]'' - T^' = D - Do (E8) 
which can be rewritten explicitely, after premultiplication by T-^^q and postmultiplication by T4, as: 

T4(0, 0; {0, 2/i}) - T4(0, 0; {0, -Eb}) = ^ ^ Ti{0, q; {0, -Eb}) [T2{q)T2{-q) - T^{q)T^{-q)] T^iq, 0; {0, -Eb}) (E9) 

q 

where in the last term we have been allowed to replace 2/i by Eb at the order we are working. The same substitution 
is allowed in Eq. (|E4[) . which completes our proof. 

Let us finally note that, with respect to this part, our procedure has obvious strong analogies with the work of 
Beliaev [23|. Indeed, as soon as we have established the non trivial result that, at our level of approximation, the 
structure of A4 can be omitted which makes it analogous to an effective interaction, our dimers behave as elementary 
bosons and we are basically back to the situation investigated by Beliaev. 
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